Linear response results for phonons and electron-phonon coupling in hep So 
fluctuations and implications for superconductivity 
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We present a FP-LMTO (full potential linear muffin-tin orbitals) study of the variation in the 
electronic structure, phonon frequencies and electron-phonon coupling in hep Sc under pressure. 
The electron-phonon coupling constant A is found to increase steadily with pressure in the hep 
phase, until the pressure reaches a value where the hep phase becomes unstable. Linear response 
calculations for the normal pressure c/a ratio predict a phase change somewhere between calculated 
pressures of 22 and 30 GPa. The calculated frequencies for the equilibrium hep lattice parameters are 
in good agreement with the inelastic neutron scattering results. There is a small upward shift in the 
F-point E2g mode frequency under pressure, in qualitative agreement with the Raman spectroscopy 
study of Olijnyk et al. (J. Phys.:Condens. Matter, 18, 10971 (2006)). From the measured value of 
the electronic specific heat constant and the calculated values of the Fermi level density of states and 
electron-phonon coupling constant, we conclude that the electron-paramagnon coupling constant in 
hep Sc should be comparable to the electron-phonon coupling constant. This indicates that the spin 
fluctuation effects are strong enough to suppress superconductivity completely in hep Sc. We argue 
that spin fluctuations should be reduced by a factor of two or more in the high pressure Sc-II phase. 
Based on estimates of the electron-paramagnon coupling constants and the calculated or estimated 
electron-phonon coupling constants, we argue that the hep phase may become superconducting with 
a very low transition temperature immediately prior to the transition to the Sc-II phase and that the 
Sc-II phase should indeed be superconducting. The electronic, electron-phonon and spin fluctuation 
properties in hep Sc under pressure are compared with those of the high pressure hep phase of Fe, 
which was reported to be superconducting a few years back. 

PACS numbers: 74.25.Jb, 74.25. Kc, 74.70.Ad, 74.90. -fn 
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I. INTRODUCTION 

This work was motivated by a recent work of Hamlin 
and Schilling,^ who reported the measurement of the su- 
perconducting transition temperature Tc in Sc as a func- 
tion of pressure. Superconductivity in Sc is induced un- 
der pressure, with T^ increasing monotonicaUy to 8.2 K at 
74.2 GPa. Hamlin and Schilling^ report measurements of 
Tc between ~ 55 GPa (T^ ~ 5K) and 74.2 GPa [T^ = 8.2 
K). The normal pressure phase of Sc is hep and it is 
known to undergo several changes in its crystal struc- 
ture as a function of pressure, the first such transition 
known to be occurring around 22-23 GPa^'— '— from hep 
(Sc-I) to a complicated structure, referred to as Sc-II. To 
date the normal and low pressure hep phase is known to 
be non-superconducting. Earlier work by Wittig et al^ 
showed that superconductivity is possibly induced in the 
hep phase at high pressure around 20 GPa, immediately 
before it enters the complex Sc-II phase. Tc in the hep 
phase was estimated to be less than 0.1 K. 

Ah initio theoretical studies of superconductivity in Sc 
as a function of pressure is rendered difficult by the fact 
that the superconducting Sc-II phase is not only com- 
plex, but also that its exact structure still remains open 
to investigation and refinementf?'^'^ Ormeci et al.^^ have 
studied the electronic structure of the Sc-II phase pro- 
posed by Fujihisa et al^ and McMahon et al.^ and pro- 
duced results suggesting that the structure proposed by 
McMahon et al. might provide a better representation 
of the Sc-II phase. Both of these proposed structures are 



composite incommensurate structures, consisting of host 
and guest substructures. 

In this communication we present ab initio calculations 
of the electronic structure, clastic properties, phonons 
and electron-phonon coupling in the low pressure hep 
phase of Sc. The purpose is to shed some light on why 
hep Sc is not superconducting, while the substance might 
become superconducting under pressure. Several scenar- 
ios are possible. It might be that the electron-phonon 
coupling in the hep phase is small, making superconduc- 
tivity unlikely in this phase. This coupling may increase 
with pressure, but still remain too small for superconduc- 
tivity to appear over the entire pressure range for which 
the hep phase is stable. Another possibility is that the 
electron-phonon coupling is strong enough to support su- 
perconductivity in the hep phase. In this latter scenario 
one needs to explore why superconductivity is suppressed 
in the hep phase and appears only in the high pressure 
Sc-II phase. Calculations of the phonon frequencies as a 
function of pressure should also show softening of certain 
modes leading to the instability of the hep phase at high 
pressure. 

We use the full potential linear muffin-tin orbital (FP- 
LMTO) method^ to study the electronic structure and 
the linear response scheme developed by Savrasov-?*- to 
compute the phonon-frcquencies, the Eliashberg spectral 
function and the electron-phonon coupling constant A as 
a function of volume in hep Sc. For the sake of simplicity 
and convenience in the calculation, the c/a ratio is kept 
fixed at the normal pressure value ji^ This restriction of 



the cja ratio being kept fixed at the normal pressure 
value should not be of any concern from the viewpoint 
of the main results and conclusions of this work. This 
will become clear from the discussion in the following 
sections. 

Sc is the lightest of the Sd- transition metals. It would 
be of interest to see how the electron-phonon coupling in 
hep Sc compares with that of the late transition metal 
hep Fe, which was reported to be superconducting by 
Shimizu ei al^ (see also RefUll) a few years back. The 
hep phase is the stable phase of Fe at pressures ^ 10 
GPa and higher. Bose et ol^ reported FP-LMTO hnear- 
response results for hep Fe. They showed, in agreement 
with earlier LAPW(linear augmented plane-wave)-based 
rigid muffin-tin (RMT) results of Mazin and coworkers;^ 
that not only should hep Fe be superconducting, but also 
conventional (electron-phonon) s-wave superconductivity 
in hep Fe should persist up to a much higher pressure 
than what is found in the experiments. Both ferromag- 
netic and antiferromagnetic spin fluctuation effects were 
considered in the works of Bose et ali^ and Mazin et 
alJ^ The conclusion was that such spin fluctuation ef- 
fects in hep Fe could lower Tc somewhat, but could not 
account for the rapid disappearance of Tc with increas- 
ing pressure, as observed by Shimizu et alJ^ Nontrivial 
differences between the electronic structures of Fe and 
Sc, due in particular to the large difference in the num- 
ber of d-electrons, cause significant differences in their 
elastic as well as electron-phonon scattering properties, 
which should have opposite effects on the overall electron- 
phonon coupling constant. 

Ormeci et al^ have used crystalline approximants 
of the incommensurate structures proposed by Fujihisa 
et a/4 and McMahon et al.^ in studying the electronic 
structure of the Sc-II phase. The unit cells of these ap- 
proximants contain 22 atoms for the model proposed by 
Fujihisa et al^ and 42 atoms for the model suggested by 
McMahon et al..— Because of the high computational de- 
mand due to the large number of atoms in the unit cells 
of these approximants, we are unable to extend the lin- 
ear response calculations of phonons and electron-phonon 
coupling to the Sc-II phase. 



II. ELECTRONIC STRUCTURE 

FP-LMTO results of total energy as a function of the 
lattice parameter in hep Sc is shown in Fig[TJ The c/a 
ratio was kept fixed at the normal pressure value. ^ As in 
the earlier calculation for hep Fcj - the generalized gra- 
dient approximation of Perdew and Wang (GGAl)ii for 
the exchange-correlation potential was used. The elec- 
tronic structure was computed using a two-K spd LMTO 
basis for the valence band. 3s- and 3p-semi-core states 
were treated as valence states in separate energy win- 
dows. The charge densities and potentials were repre- 
sented by spherical harmonics with I < 6 inside the non- 
overlapping MT spheres and by plane waves with ener- 
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FIG. 1: FP-LMTO total energy per atom versus lattice pa- 
rameter in hep Sc for the c/a — 1.592, valid for the normal 
pressure values: a — 3.31 A, c = 5.27 A[y]. For convenience in 
plotting, a constant value of -1529.0 Ry from the total energy 
per atom has been subtracted. 



gies < 48-70 Ry, depending on the lattice parameter, in 
the interstitial region. Brillouin zone (BZ) integrations 
were performed with the full-cell tetrahedron method^ 
using 1200 k-points in the irreducible zone. The pressure 
and bulk modulus calculated from this energy-volume 
curve, using the generalized Birch-Murnaghan equation 
of state^^iii are shown in Table HI Note that the mini- 
mum energy occurs almost exactly at the normal pressure 
lattice parameter of 6.255 a.u. The calculated bulk mod- 
ulus of 56 GPa at this^^ equilibrium lattice parameter is 
somewhat higher than the experimental value, 43.5 GPa. 
One interesting observation is that the lattice parameter 
5.75 a.u., or the volume per atom Vq = 131 a.u., corre- 
sponds to a pressure of 22 GPa. This is the pressure at 
which the hep Sc-I phase becomes unstable against the 
formation of a complicated structure known as the Sc-II 
phasei^i^i^ As will be discussed later, the linear response 
results show that the phonon frequencies calculated for 
lattice parameter a — 5.6 a.u. are imaginary. Thus, for 
the normal pressure c/a ratio, the hep structure becomes 
unstable at some lattice parameter between 5.75 and 5.6 
a.u. The calculated pressure at 5.6 a.u. is ^ 30 GPa. 
Given the uncertainties in the pressure calculation, this 
result is in good agreement with the experimental obser- 
vation of a phase change in Sc at 22-23 GPa. Additional 
linear response calculations to better locate the phase 
change were not done. 

The FP-LMTO DOS calculated for three lattice pa- 
rameters are shown in Figl^l Panel (a) shows the DOS 
for the experimental normal pressure lattice parameters: 
a = 3.31 A, c = 5.27 A,^ Panels (b) and (c) show 
the DOSs for the lattice parameters 5.75 a.u. and 5.60 
a.u., respectively, with the same c/a ratio as in (a). For 
the equilibrium lattice parameters (panel (a)), the Fermi 
level falls in a shallow valley between two small peaks. 



TABLE I: FP-LMTO results for hep Se for the normal pres- 
sure cja ratio, 1.592. &= lattice parameter (a.u.), Vo = vol- 
ume per atom (a.u.), P=pressure (GPa), B^Bulk Modulus 
(GPa), A''(0)= DOS at the Fermi level (states/(Ry cell)). 



TABLE II: LMTO-ASA results for (approximate) s-, p-, d- 
and /-orbital resolved charges: ns,np,nd,nf and the corre- 



a 
Vo 
P 
B 



6.255 6.15 

168.72 160.37 

-0.084 3.21 

56.2 67.2 



6.05 

152.67 

6.85 

79.6 



5.95 

145.22 

11.22 

93.0 



5.85 5.75 

138.02 131.07 

16.20 22.28 

108 125 
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FIG. 2: FP-LMTO DOSs for three lattice parameters of hep 
Sc: (a) Normal pressure values[3|: a = 6.255 a.u., c/a — 
1.592, (b) a — 5.75 a.u. and (c) a = 5.60a.w.; c/a ratio for (b) 
and (c) is the same as in (a) . It appears that as the Fermi level 
moves from the shallow valley region between two small peaks 
(a) to a shoulder (c), the hep structure becomes unstable. 



With increasing pressure, the bands broaden, flattening 
these peaks and the valley region (panel (b)). It appears 
that as this region around the Fermi level changes to a 
shoulder (panel (c)), the hep phase becomes unstable. 

Hamlin and Schilling^ discuss the importance of the 
role of s — > d charge transfer in the variation of electron- 
phonon coupling and the changes in the crystal structure 
under pressure. Pettifor— has discussed the variation of 
the d-band occupancy under pressure of Ad- transition 
metals, showing that the increase in the d-electron oc- 
cupancy of the early transition metals under pressure is 
rapid. Duthie and Pettifor— used a s — d band model 
to discuss the correlation between the d-band occupancy 
and the crystal structures in rare earths. They showed 
that d-electron concentration nd increases under pressure 
and a sequence of crystal structure changes takes place 
as Ud changes from 1.5 to 2.5. With increasing pres- 
sure, the volume available to the electrons outside the 
ion cores diminishes rapidly. As a result, simple metals 
like Li and Na are known to become non-free-electron- 
like under pressure and Cs is known to become a tran- 
sition metal as the 5d-band begins to fill due to s — > d 



sponding Fermi level DOSs: Ns,Np,Na,Nj: 
n^ = ijid + nf) — {us + np) and A*'' 
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j/,j^^ are presented to 
reveal the trends in the inter-orbital charge transfer and the 
redistribution of the orbital-resolved DOSs as a function of 
volume per atom. N{0) is the total DOS at the Fermi level. 
All DOS are in units of states/ (Ry cell). &— lattice parameter 
(a.u.) 



a 

ris 

Up 

rid 

Uf 

n~ 
n' 



6.255 
0.730 
0.690 
1.545 
0.031 
0.156 
1.11 



6.15 
0.720 
0.677 
1.567 
0.033 
0.203 
1.15 



6.05 
0.714 
0.662 
1.589 
0.034 
0.247 

1.18 



5.95 
0.704 
0.645 
1.614 
0.036 
0.301 

1.22 



5.85 
0.695 
0.626 
1.641 
0.037 
0.357 

1.27 



5.75 
0.685 
0.603 
1.672 
0.039 
0.423 

1.33 



A(0) 
A,(0) 
A,(0) 
A^d(O) 

N' 



57.89 
0.941 
14.836 
40.789 
1.326 
2.67 



54.78 
0.946 
13.729 
39.110 
1.301 
2.75 



53.03 
0.964 
12.889 
37.887 
1.287 
2.83 



51.60 
0.982 
12.257 
37.080 
1.280 
2.90 



50.16 
1.020 
11.914 
36.444 
1.280 
2.92 



53.34 
1.226 
13.558 
37.222 
1.335 
2.60 



charge transfer (see references in [1|). In any electronic 
structure calculation the numbers of s-, p- and d- elec- 
trons are essentially functions of the basis set used, if 
such a division is at all possible for the method used. 
It is easier to keep track of s — > d charge transfer in 
LMTO-ASA (atomic sphere approximation),^^- where the 
basis consists of muffin-tin orbitals only, rather than in 
FP-LMTO. In Table [U we present the orbital- or partial 
wave channel-resolved electron numbers as a function of 
the lattice parameter, as obtained by the LMTO-ASA 
methodi^ The results presented were obtained by using 
the exchange correlation potential of Perdew and Wang— 
in the local density approximation. Checks for a couple 
of lattice parameters using GGAlii had revealed similar 
results. The decrease in Us and the increase in Ud are 
monotonic as a function of decreasing lattice parameter. 
In fact, charge transfer takes place from the extended s- 
and p-orbitals to the less extended or localized d- and 
/-orbitals. Of particular interest are n' and n~, defined 

'"^"^ '^' = n!+n^ ^"^"^ "" = ("-<* + "■/) ~ ('^s + "p)' quantify- 
ing the transfer of electrons from the delocalized to the 
localized channels. At normal pressure the electrons are 
divided almost equally between the localized and delocal- 
ized channel (n~ = 0.15). With increasing pressure the 
electrons more and more occupy the localized orbitals. 

As the pressure increases different partial bands 
broaden at different rates. In Table [Til we also show the 
changes in partial or angular momentum-resolved DOSs 
at the Fermi level. The ratio of the localized orbital 
DOSs (iVd(O)-t- Nf{Q)) to the delocalized orbital DOSs 
{Ns{Q)+Np{Q)) seem to increase steadily under pressure. 



except for the lowest lattice parameter shown. The de- 
crease m N' for the lattice parameter 5.75 a.u. could be 
specific to LMTO-ASA, as the somewhat large change in 
A^(0) between the lattice parameters 5.85 and 5.75 a.u. 
is not observed in the FP-LMTO results (Table |T|. 

Ormeci et alM- have calculated the partial charges in 
the high pressure Sc-II phase. They use the full-potential 
local-orbital (FPLO) method of Koepernik and Eschrig^'* 
with an spd-basis. At ambient pressure, in the hep struc- 
ture, their nZ-projected charges for the 4s, Ap, and 3d 
channels are 0.70, 0.59 and 1.70, respectively. In the 
Sc-II phase at about two thirds of the ambient pressure 
volume per atom, these numbers become 0.54-0.65, 0.27- 
0.51, and 1.92-2.05, respectively. Their calculations also 
suggest that a complete s ~> d charge transfer would 
require pressures in excess of 240 GPa. 



III. PHONONS AND ELECTRON-PHONON 
COUPLING 



The FP-LMTO linear response results for phonons and 
electron-phonon coupling are summarized in Table IIIII 
The dynamical matrix was generated for 32 phonon wave 
vectors in the irreducible BZ, corresponding to a mesh 
of (7,7,7) reciprocal lattice divisions. The BZ sums for 
the dynamical and electron-phonon (Hopficld) matrices 
was done for a (28,28,28) mesh, resulting in 1200 wave 
vectors in the IBZ. These choices were based on exten- 
sive tests performed by Savrasov and Savrasov^'* in their 
study of the phonons and electron-phonon coupling in 
elemental metals and alloys. Phonon density of states 
F{ui), Ehashberg spectral function a'^F{u!) and the func- 
tion a^(w) = a^F(uj)/F{uj) are shown in Fig. [3]for three 
different lattice parameters with the c/a ratio fixed at 
the normal pressure value. 

The calculated phonon frequencies appear to be in rea- 
sonably good agreement with the inelastic neutron scat- 
tering (INS) results.-^- "^^ For the normal pressure lattice 
parameter (6.255 a.u.), the calculated longitudinal and 
transverse optic (LO and TO) mode frequencies at the 
F-point are 6.273 THz (209 cm"!) and 3.912 THz (130.5 
cm~^), respectively. These compare well with the inelas- 
tic neutron scattering results of Wakabayashi et alM: 
6.91 and 4.04 THz for the LO and TO modes, respec- 
tively. The recent Raman spectroscopy study of Olijnyk 
et alM- puts the TO mode {E2g) frequency at a somewhat 
higher value of 139 cm~^. On the whole, the calculated 
frequencies are a bit lower than the INS results of Wak- 
abayashi et alM These authors have analyzed the INS 
results in terms of a sixth-neighbor modified axially sym- 
metric force constant model and produce a frequency dis- 
tribution, which reveals a maximum frequency of about 
7.25 THz or 241.8 cm~^, approximately 10% higher than 
the maximum calculated frequency 215 cm~^ f Table HIT)) 
in this work. The shape of the frequency distribution 
agrees in general with the phonon DOS F(uj) shown in 
Fig. [31 except for the relative heights of some of the low 



frequency peaks. 

The Raman spectroscopy work by Olijnyk et al^ 
shows a shift of the E2g frequency under pressure from 
139 cm^i to about 150 cm"! at 18.8 GPa. The cal- 
culated F-point TO mode frequency increases from 130.5 
cm^^ at the normal pressure lattice parameter (a= 6.255 
a.u.) to 135.1 cm~^ at a = 5.75 a.u., i.e., the increase 
is much more subdued. In Table IIIII we have presented 
the calculated E2g frequency for various lattice parame- 
ters. Although the variation may seem non-monotonic, 
the differences are very much within the error bars of the 
calculation. What can perhaps be concluded is that there 
is no phonon-softening and most probably this frequency 
does increase by approximately 5 cm~^. This is about 
45-50% of the observed shift. Note that calculations for 
the stable hep phase for a lattice parameter between 5.75 
and 5.6 a.u. might reveal a higher shift. The discrepancy 
between the calculated and the observed shift could be 
partly due to the possible variation in the c/a ratio under 
pressure, which has not been considered in the calcula- 
tion. Olijnyk et al^ have analyzed the shift in the E2g 
frequency in terms of the relation 
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1 M^/3a2C44 
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(1) 



where m is the atomic mass and C44 is the elastic shear 
modulus. They relate their result to the behavior of the 
shear modulus at the transition. The possible variation 
of the c/a ratio with pressure should be factored in in 
this discussion. The coarse wave vector mesh used for 
the phonon frequency calculation prevents us from us- 
ing the slope of the phonon dispersion curves or the long 
wavelength method to compute the elastic constants. Ab 
initio calculation of the elastic constants, including C44, 
via long wavelength and homogeneous deformation meth- 
ods using energy-minimized c/a ratios will be the subject 
of a separate publication. 

It is noteworthy that both calculation and the Raman 
study show a positive shift of the E2g frequency, in con- 
trast to the behavior in Y and the regular Lanthanides, 
where a softening of the i?2g mode under pressure is asso- 
ciated with the transition to the high-pressure Sm-type 
structurc i^^'^^ Incidentally, in our linear response calcula- 
tion, as the lattice parameter is lowered to 5.6 a.u. where 
the instability of the hep phase is manifest via appearance 
of some imaginary frequencies, the E2g mode frequency 
stays real and has a lower value 124 cm~^. At a lattice 
parameter of 5.2 a.u., this frequency, along with many 
others, are imaginary. 

Broadening of the phonon bands and the increase in 
the maximum phonon frequencies with pressure is appar- 
ent in Fig. [3l Maximum phonon frequency increases from 
215 cm^^ to 295 cm~^, as the lattice parameter changes 
from the normal pressure value to 5.75 a.u., where the 
estimated pressure (Table |l| is 22-23 GPa. Note that 
with increasing pressure, initial increases in w, Co and wi„ 
are followed by a decrease, revealing phonon softening. 
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FIG. 3: Phonon density of states F{lo), the Eliashberg spectral function a F{uj) and the function a (uj), defined as the ratio 
o?F{ui)/ F{u)), for three lattice parameters a — 6.255 (experimental normal pressure value), 6.05, and 5.75 a.u., with c/a ratio 
fixed at the normal pressure value. 



At normal pressure, the shape of the a^F{uj) function 
follows that of the phonon DOS F{lo). With increas- 
ing pressure, differences between the shape of the two 
functions appear due to increased contribution from low 
frequency phonons, a consequence of phonon softening. 
There is a feed-back effect: electron-phonon coupling 
leads to both phonon linewidth and rcnormalization of 
the phonon frequency, and the latter, in turn, affects the 
coupling constant. Unfortunately, the 7,7,7 division of 
the BZ misses most of the symmetry points. However, it 
appears that there are several regions of the wave vector 
space that show large mode coupling. Linear response 
calculations for the lattice parameter 5.6 a.u. show that 
the phonon frequency at the wave vector closest to the 
L-point for the 7,7,7 division becomes imaginary first. 
For the lattice parameter 5.75 a.u., this wave vector has 
a very large mode coupling constant. In Table Hill it is 
clear that general phonon softening starts as the lattice 
parameter increases beyond 6.05 a.u. (as seen from w;„) 
or 5.85 a.u. (as seen from ui and w), while the £'25 mode 
continues to stiffen and the maximum phonon frequency 



continues to increase. 

The electron-phonon coupling parameter \ph is a com- 
bination of an electronic parameter (Hopfield parame- 
ter) ry = (/^)A^(O) (A^(0) being the Fermi level DOS 
for one type of spin) and the mean square phonon fre- 
quency (oj^): Xph — rj/m{Ld'^). (/^) is the Fermi surface 
average of the square of the electron-phonon matrix el- 
ement. The electronic and phonon-related parameters 
act in opposite directions in affecting the coupling con- 
stant: A is enhanced by having higher {P) at lower fre- 
quency. As the Hopfield parameter for transition metals 
has often been calculated using the rigid muffin-tin ap- 
proximation of Gaspari and Gyorffy,'^^ in Table IIIII we 
have compared the values obtained via the FP-LMTO 
linear response method (77) and those obtained by using 
the RMT scheme implemented within the LMTO-ASA 
method, -^-i^ known as the rigid atomic sphere (RAS) 
method. The latter values, labeled as ri{RMT / RAS) , are 
in general agreement with the FP-LMTO results, but are 
somewhat overestimated, with the level of the overesti- 
mation increasing with pressure. This is in contrast with 



the results for the late transition metal hep Fe, where 
RAS was shown to consistently underestimate the Hop- 
field parameter, with the level of the underestimation in- 
creasing with pressure.— The values of the Fermi surface 
averaged electron-phonon matrix element (/^), obtained 
from the linear response calculation and by using the 
RMT approximation, are also shown in Table Hill 



IV. SUPERCONDUCTIVITY IN HCP Sc, 

BASED PURELY ON ELECTRON-PHONON 

COUPLING 

Since ab initio results of the electron-phonon coupling 
in hep Sc have not been reported in the literature, it 
would be of some interest to derive values of the su- 
perconducting transition temperature Tc based on the 
present calculations. Indeed, the values of the coupling 
constant Xph listed in Table IIIII suggest that Tc could 
be high enough to be experimentally observable. The 
superconducting transition temperature can be obtained 
by solving the linearized isotropic Eliashberg equation at 
Tc(see, e.g., ReiM): 

Z{iio^) = l + ^yW^+(n-n')signK0(2) 



n' 



Z{iiu„)A{iujn) = ttTc J2 ^-(^ 



, A{iUn') 



wherein — 7rTc(2n+l) is a Matsubara frequency, A(ia;„) 
is an order parameter and Z{iLUn) is a renormalization 
factor. Interactions W+ and W- contain a phonon con- 
tribution Xph, a contribution from spin fluctuations A^/, 
and effects of scattering from impurities. With scattering 
rates 7™ = t^ ^^'^ Inm = -2^ — referring to magnetic 



and nonmagnetic impurities, respectively, the expressions 
for the interaction terms are : 

W+{n — n') — Xph{n — n') + Xsf{n — n') 

+Snn'il7im+lm), (3) 



and 



W-{n — n') = Xph{n — n') — Xsf{n — n') 

-fl*{uJc) + 6nn' ilnm -7m)- (4) 



The phonon contribution is given by 

dujuja^{uj)F{uj) 



Xph{n - n') = 2 



(W„ -OJn'Y +W^ 



(5) 



where a^(w)F(w) is the Eliashberg spectral function, de- 
fined as 

^ ' k,k',y> 

(6) 



Here, g^ k' i^ ^^^ electron-phonon matrix element, with 
u being the phonon polarization index and k, k' repre- 
senting electron wave vectors with band indices i, and j, 
respectively. Ap;i(0) = Xph is the electron-phonon cou- 
pling parameter, the values of which are given in Table 
mil The contribution connected with spin fluctuation can 
be written as 



Xsf{n - n') 



duj^P{uj) 



{uJn-^n'Y +^^ 



(7) 



where P{io) is the spectral function of spin fluctuations, 
related to the imaginary part of the transversal spin sus- 
ceptibility x±i^) ^^ 

F(tj) =--(|5kk'|^Imx±(k,k',cj)) 

TT \ I FS 

where ( )ps denotes Fermi surface average. Xsf — As/(0) 
is often referred to as the electron-paramagnon coupling 
constant. 

In Eq.(j4]), /i*(ijJc) is the screened Coulomb interaction. 



^i*{uJc) = 



M 



l + ^ln(£;/wc)' 



(8) 



with 1.1 = {N{Q)Vc)pg being the Fermi surface average 
of the Coulomb interaction. E is a characteristic electron 
energy, usually chosen as the Fermi energy Ep and ujc is a 
cut-off frequency, usually chosen ten times the maximum 
phonon frequency: ujc — lOw^Jf'^. 

For a start, we ignore all consideration of spin fluc- 
tuations and impurity scattering and solve the Eliash- 
berg equation with only the electron-phonon term and 
the Coulomb pseudopotential ij,*{uJc)- As is often done, 
we assume that a reasonable value for ^ is ~ 1.0, and 
from the calculated Fermi energies Ep we obtain /x* for 
all volumes, with the cut-off frequency uJc assumed to be 
ten times the maximum phonon frequency. The values 
are listed in Table Hill As shown in this table, the super- 
conducting transition temperature in hep Sc based on 
consideration of electron-phonon coupling alone can be 
significant, increasing monotonically with pressure from 
2K to ~ lOK until the instability of the hep phase sets 
in. Interestingly, the Tc values reported by Hamlin and 
Schilling^ for the high pressure Sc-II phase fall in this 
range, while no superconducting behavior has been ob- 
served in the hep phase. 

For pedagogical reasons, we have listed in Tabl dllll the 
values of Tc obtained by using the Allen-Dynes form"^" of 
the McMillan expression: 



Tc^ 



1.2 



exp 



1.04 (1 + Ap,,) 
A-/i*(l-f 0.62Apft) 



(9) 



where coin is the logarithmically averaged phonon 
frequency;^ obtained from our linear response calcula- 
tions and reported in Table HTTl Note that the Coulomb 
pseudopotential /i* appearing in the McMillan equation 



TABLE III: Hopfield parameters from the linear response calculation rj and the rigid muffin-tin (atomic sphere) approximation rj 
(RMT/RAS), the Fermi surface average of the square of the electron-phonon matrix element (I^) and (/^)(RMT/RAS) from the 
linear response and the RMT results, respectively; logarithmic average phonon frequencies coin; maximum phonon frequencies 
LUm, average phonon frequencies to — {co) and iD = (oj^)^'^; the F-point i?2g mode frequency u{E2g)\ Coulomb pseudopotentials 
for Eliashberg equation (/i*(cJc)) and McMillan formula (/i*n); electron-phonon coupling parameters \ph', calculated critical 
temperatures (T^''''^) from the solution of the Eliashberg equations ((2)1 and the critical temperatures from the McMillan formula 



®(r,^ 



for various lattice parameters a. 



a 


Ob 


6.255 


6.15 


6.05 


5.85 


5.75 


n 


Ry/bohr^ 


0.0381 


0.0433 


0.0519 


0.0666 


0.0771 


{!') 


(Ry/bohr)2 


0.0026 


0.0031 


0.0040 


0.0056 


0.0065 


?7 (RMT/RAS) 


Ry/bohr^ 


0.0439 


0.0509 


0.0591 


0.0826 


0.1077 


(/2) (RMT/RAS) 


(Ry/bohr)2 


0.0030 


0.0037 


0.0045 


0.0065 


0.0081 




K 


160.4 


195.0 


198.1 


195.4 


177.7 


^In 


cm^^ 


111.5 


135.5 


137.7 


135.8 


123.5 


(^m 


cm~^ 


214.5 


230.4 


246.9 


279.6 


294.7 


U[E2,) 


cm~ 


130.5 


134.95 


133.4 


136.14 


135.0 


LU 


cm~ 


122.4 


142.3 


145.1 


145.6 


135.3 


LU 


cm~^ 


131.5 


148.6 


152.3 


155.8 


148.1 


M*(^c) 




0.253 


0.256 


0.260 


0.266 


0.269 


^i*{uJlu ) 




0.160 


0.161 


0.163 


0.165 


0.166 


\h 




0.639 


0.576 


0.657 


0.807 


1.033 


j.McM 


K 


2.34 


1.66 


2.96 


5.66 


9.18 


T-Tcalc 


K 


2.17 


1.75 


3.13 


6.00 


9.64 



above is related to ^*{ujc) appearing in the Eliashberg 
equation via— 



IJL = PL [LOin) 



lJL*{uJc) 



[1+ ^l*{uJc)^u{uJc/l^Om)) 



(10) 



Our results are computed with Lod^m = 10. The T^. 
values obtained by solving the Eliashberg equations and 
those from the McMillan expression Eq. ^ show excellent 
agreement. Earlier, calculations for hep Fei^ had shown 
the McMillan expression to overestimate T^ with respect 
to the results from the Eliashberg equation, while for fee 
and bet boron an opposite trend was revealedi^ 



V. SPIN FLUCTUATION EFFECTS 

A. hep So at normal pressure 

Faced with the results of the previous section, one 
needs to produce convincing arguments as to why hep 
Sc is not superconducting despite sufficiently strong 
electron-phonon coupling. A mechanism that is known 
to cause suppression of superconductivity is spin fluctu- 
ations, which is often quoted as the reason why fee Pd is 
not supcrconductingf2I It is argued that the high value of 
A^(0) in fee Pd leads to considerable Stoner enhancement 
of paramagnetic spin susceptibility, making it a border- 
line ferromagnetic material, and spin fluctuations tending 
to ferromagnetic alignment of spins lead to the breaking 
of the Cooper pairs. That such a mechanism is operative 



in hep Sc as well is highly probable and has been dis- 
cussed in the literature on a few occasions i^^*^ Jensen 
and Maita^ argue that the spin fluctuation effects are 
responsible for the rapid depression of T^ in the Zr-Sc 
alloy system, as Sc is added. Ah initio calculations of the 
spin susceptibility of hep Sc at equilibrium volume based 
on the spin-density functional theory by MacDonald et 
al^ yield an exchange-correlation enhancement factor of 
4.03 over the band value and 17.2 over the free-electron 
value at the same average electron density as in hep Sc. 
An earlier calculation by Das^ puts the Stoner enhance- 
ment factor at 4.6 (over the band value). The value cal- 
culated for fee Pd by Janak'*- is 4.46. This shows that 
the spin fluctuations in hep Sc should be as strong as in 
fee Pd. 

For a proper theoretical treatment of the spin fluctu- 
ation effects one needs to compute Xsf{n — n') from the 
spin susceptibility function given by Eq.(I7]). However, it 
is important to note that such treatments tacitly assume 
a Migdal-like theorem being applicable to spin fluctua- 
tions. The Eliashberg equations (Eq.©) are based on 
the assumption that the maximum or the cut-off energy 
of spin fluctuation is much smaller than the characteristic 
electronic energy, e.g. the Fermi level. A somewhat qual- 
itative treatment of spin fluctuations can be based on es- 
timating \sf — As/(0) from experiments. Both electron- 
phonon and the electron-paramagnon interactions con- 
tribute to the electronic specific heat. In an independent 
one-electron picture this is interpreted as the electronic 
mass enhancement or equivalently, enhancement of the 
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density of states over the bare value N{Q). The latter is 
the value given by calculations, where these interactions 
are not included in the one-electron Hamiltonian. Thus, 
a reliable estimate of the electron-paraniagnon coupling 
constant As/ can be obtained from the measured value of 
the temperature co-efRcient of the electronic specific heat 
7, and the calculated values of the bare electron density 
of states and the electron-phonon coupling constant Xph'- 



1 = 


= y4^*(0) 


(11) 


iV*(0) = 


= iV(0)(l + Ae//) 


(12) 


^eff = 


= Ap/i + Xgf ■ 


(13) 



Here, 7 and N{0) refer to the values per atom. The 
Coulomb interactions are included in an average sense 
in the density functional calculations of N{0), and have 
therefore been left out of Eq. (12). 

Among all the elemental metals, excluding the rare- 
earths, Sc has the largest electronic specific heat con- 
stant 7, followed by Y and Pdj^^ The reported experi- 
mental values at normal pressure are 10.9-10.33 mJ/mole 
j^2 ^^^42^44 Considering the latest and the most 
conservative value of 7 =10.33 mJ/molc K^,^^'^^ we get 
^'*(0)= 2.2 states/(eV atom spin). With the calcu- 
lated value of 1.067 states/ (eV atom spin), we get Ae// 
= 1.063. The calculated value of Xph — 0.639 then 
yields Xgf =0.422. Effects of spin fluctuations on Tc 
can be incorporated by a simple rescaling of Xph and 
the Coulomb pseudopotential /i*: Xph -^ Xph /{I + A^/), 
/i* ^ {fi* + Xsf)/{l + Xsf).'^^ An extension of the McMil- 
lan formulaiS that is often used to incorporate the spin 
fluctuation effects is 



ph 



Tr 



1.2 



exp 



1.04(1 + Ap,, + X,f) 



Xph - Xsf - Ai*[l + 0-62{Xph + As/)] 



This formula is meaningful as long as As/ is sufficiently 
less than Xph, so that the denominator in the argument 
of the exponential in Eq. (fTi]) stays positive and not 
close to zero. For the above values of Xph, As/ and the 
Coulomb pseudopotential /z*=0.16 (see Table HIT)) , this 
condition breaks down. In fact, for comparable values 
of Xph and As/, the Kernel of the Eliashberg equation 
(Eq. (21)) W-{n — n'), given by Eq.Q, becomes negative, 
allowing no solution for Tc- It should be emphasized 
that the McMillan-type formulas in the presence of spin 
fluctuations are good for rough ball-park estimates of Tc 
only. 

Note that our calculated value of Xph = 0.639 for nor- 
mal pressure hep Sc is significantly higher than the value 
Xph = 0.3 that has been suggested by Knapp and Jones, " 
based on a comparison of the high and low temperature 
electronic specific heats. A value of Xph = 0.3 would sug- 
gest Xsf — 0.76, more than twice larger than Xph- Our 
results suggest that As/ < Xph, but their values are close. 
This leads to the possibility that if spin fluctuations can 
be suppressed via application of pressure, superconduc- 
tivity can indeed appear. 



Based on the analysis of high and low temperature 
specific heats, Knapp and Jones^ also suggested that 
for Pd Xph — 0.7. The electronic specific heat constant 
7 for Pd is only slightly lower than that in hep Sc^^i^ 
lying between 10.0 and 9.2 mJ/mole K^. The calculated 
bare band density of states A^(0) for fee Pd is about 1.3 
states/(eV atom spin), giving Ae// in the range 0.66 and 
0.5. An analysis by Savrasov and Savrasov§ puts the 
value at 0.69. The value Xph = 0.7 suggested by Knapp 
and Jones would leave no room for spin fluctuation effects 
in Pd, and would render Pd superconducting. The FP- 
LMTO linear response calculation for fee Pd by Savrasov 
and Savrasov^ yields Ap^ = 0.35, giving As/ in the range 
0.3-0.15. These values of Xph and As/ can adequately 
explain the nonexistence of superconductivity in fee Pd. 
We have repeated the linear response calculation for fee 
Pd and obtain a value of Xph = 0.37, in close agreement 
with the result of Savrasov and Savrasov.® Note that the 
calculation by Pinski and Butler^^ based on the RMT ap- 
proximation yields Xph = 0.41 for fee Pd, in reasonable 
agreement with the FP-LMTO linear response results. 
To summarize, our linear response calculations and elec- 
tronic specific heat analysis show that in both hep Sc 
and fee Pd the values of electron-phonon and electron- 
paramagnon coupling constants are comparable, with the 
latter being slightly lower. This is qualitatively different 
from the previous results of Knapp and Jones ,^^ which 
would suggest that in Sc Xph is significantly lower than 
As/, and vice versa for fee Pd. Note that according to 
our results both Ap^ and As/ are larger in hep Sc than 
the corresponding quantities in fee Pd. 

While the nature of spin fluctuations in fee Pd is ac- 
cepted to be ferromagnetic, that in hep Sc could be an- 
tiferromagnetic. The dynamic susceptibility x((7, lo) in 
Pd has a peak at g = 0, whereas xil^'-^) in hep Sc is 
expected to have a peak at some finite wave vector Q. 
First-principles calculation of xil) including spin-orbit 
coupling and all other relativistic effects by Thakor et 
alM- yields a peak at the wave vector q = (0, 0.0.57)7r/c), 
consistent with a Fermi surface nesting vector. Earlier 
calculations produce a peak along the same direction 
but at different values^i^ of q. Capellmann^ has dis- 
cussed the effect of incipient antiferromagnetism in Sc 
and shown that the antiferromagnetic spin fluctuations 
should lead to a repulsive electron-electron interaction, 
resulting in suppression of superconductivity. The mod- 
ified McMillan formula (Eg. p^ 'l has been used to es- 
timate Tc for both ferromagnetic and antiferromagnetic 
spin fluctuationSfi^ii^ 



B. Superconductivity in hep Sc at high pressure 
immediately before the transition to Sc-II phase 

Our linear response calculation yields a high value 
Xph ^ 1.0 of the electron-phonon coupling constant at 
high pressure hep phase prior to the transition to the 
complex Sc-II structure. It is natural to ask whether 



spin fluctuations should be able to suppress supercon- 
ductivity despite such strong electron-phonon coupling. 
Since specific heat data are not available for such high 
pressures, a reliable estimate of Xgf is difficult. However, 
a reasonable step might be to scale the ambient pres- 
sure As/ according to the density of states N{0) at high 
pressure shown in TableU This gives Xgf = 0.346 for the 
smallest volume or the highest pressure calculation. With 
Xph = 1.0, As/ = 0.346, /I* = 0.166, and ujin = 177.7X, 
as shown in Table lllli and using Eq. p^ , we get T^ — 
0.14 K. As pointed out earlier, at comparable values of 
Xph and As/, Eq. p^ may not be reliable, and one needs 
to solve the Eliashberg equation. However, our estimate 
of Tc = .lAK indicates that there is a possibility of hep 
Sc turning superconducting just before the onset of the 
high pressure Sc-II phase. This is indeed what has been 
reported by Wittig et al.^ Note that such a conclusion 
would not be tenable for smaller values of Xph, as the 
results of Knapp and Joncs*^ would suggest. 

The use of energy-optimized c/a ratio at higher than 
ambient pressures may change the calculated values of 
Xph somewhat, but such small changes are not expected 
to influence our results qualitatively and will not affect 
the nature of the conclusions. 



C. Superconductivity in Sc-II phase 

Spin fluctuations are expected to be significantly re- 
duced in the Sc-II phase. According to the calculations 
of Ormeci et al.^ the density of states at the Fermi 
level N[Q) in the Sc-II structure should be about half the 
value in the ambient pressure hep Sc. This suggests that 
As/ should be reduced to a value of 0.2 or less. Reason- 
able choices of lj;„ in the Sc-II phase at pressures ~30 
GPa (see TableHni and ^l* should be 210.0 K and 0.166, 
respectively. A value of Ap/i=0.95 in Eq. p^ would yield 
a rc=l-4 K. Reducing the value of As/ to 0.1 would result 
in a Tc of 4.5 K. These numbers suggest that supercon- 
ductivity in the Sc-II phase can be explained based on 
the quantities calculated for the high pressure hep Sc and 
reasonable estimates of the reduction in the spin fluctu- 
ation effects that is expected as a result of the change in 
the density of states. Note that phonon softening causes 
Xph to increase in the hep phase immediately before the 
transition to the Sc-II structure. At the start of the Sc-II 
phase the value of Xph should be somewhat less than the 
value 1.033 shown in Table [1111 hence the choice 0.95. It 
should be noted that although the density of states A^(0) 
decreases with increasing pressure, there is a rapid in- 
crease in the electron-phonon matrix element (/^) with 
decreasing volume (see Table IV). A high value of Xph 
in the Sc-II phase is thus possible despite the reduced 
value of -/V(O). Chances are that at the start of the Sc- 
II phase both Xph and As/ are somewhat less than the 
values suggested here. 



VI. COMPARISON BETWEEN HCP Sc AND 
HCP Fe 

The bulk moduli of all transition metals are known to 
increase (at least initially) as a function of band fiUingf^ 
and the bulk moduli of the late transition metals are in 
general higher than those of the early ones. The bulk 
modulus of Fe is thus expected to be larger than that 
of Sc, a consequence of Fe having more d-electrons. The 
calculations of Bose et alM^ for the hep phase of Fe show 
values in the range of 300 GPa to 970 GPa, corresponding 
to volumes per atom for which the calculated Tc is above 
zero (the theoretical superconducting phase). These val- 
ues are 4-5 times higher than the bulk modulus values 
of Sc shown in Table HI Accordingly, the average phonon 
frequency in hep Fe is about four times higher than in 
Sc. Because of a much larger number of d-electrons in 
hep Fe (6.5 on average) than Sc (1.6 on average), the 
electron-phonon matrix element (/^) in Fe is also much 
larger. As a result, despite the lower value of iV(0), the 
Hopfield parameter 77 = N{0){P) in hep Fe turns out 
to be more than double that in Sc (compare Table IIIII 
with TABLE II of Ref. [12]). However, the dominating 
effect turns out to be the lower phonon frequencies in Sc, 
leading to higher values of Xph ■ 

Spin fluctuations, most probably of an antiferromag- 
netic natureji^ii^ are believed to be present in the high 
pressure hep phase of Fe. Because the density of states 
A^(0) in hep Fe is about half that in hep Sc, spin fluc- 
tuations in hep Fe should be much weaker. There is no 
specific heat data available for hep Fe, the stable phase 
of Fe at pressures of ~10 GPa and higher. As a result, 
an experimental estimate of As/ in hep Fe is not avail- 
able. Some theoretical estimates were provided in Refs. 

nil. 



For a crude estimate of the ratio between the values 
of Xph for hep Fe and Sc, we can assume the Hopfield 
parameter 77 to be simply proportional to the fractional d- 
DOS, Nd{Q)/N{Q). For an estimate of m{uP'), we can use 
the result proposed by Moruzzi, Janak and Schwarz^ 
relating the Debye temperature Qd to the bulk modulus, 
atomic mass and the average Wigner-Seitz radius: 



e. 



41.63 



soB 



(15) 



where sq is the average Wigner-Seitz radius in a.u. and 
B is the bulk modulus in kbar. This result suggests that 
mjuj'^) (X sqB. Using B values from Tableland TABLE 
I of [l2|, the DOS values from Table [III sq computed 
from a — 4.6 a.u, c/a = y/sjs for Fe and a — 6.255 
a.u., c/a = 1.592 for Sc, we get Xpe/Xsc = 0.45. This 
compares favorably with the value 0.66, according to the 
computed linear response values (Table |TTT] and TABLE 
H of [12). We have used Nd{0) = 34.55 states/(Ry cell) 
and A^(0) = 40.79 states/(Ry cell) (TABLE I of [13) 
for Fe. The proportionality of 77 to fractional d-DOS, 
Nd{0)/N{0), can be somewhat justified on the basis of 
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the RMT result of Gaspari and GyorfFy^: 

, = 2iV(0)E(/ + l)M5..^|^, (16) 



where N{0) is the Fermi level DOS per atom per spin 
and fi is a relative partial state density, 



fi 



Ni jO) 
N{0) 



(17) 



M/^;+i is the electron-phonon matrix element obtained 
from the gradient of the potential and the radial solu- 
tions Ri and Ri+i of the Schrodinger equation evaluated 
at the Fermi energy. Neglecting the matrix elements and 
all partial DOSs other than Nd{Q) amounts to the result 
rj ex Nd{0)/N{0). Of course, this provides a very crude 
estimate and can only be relied on in deciding whether A 
for one material is greater or lower than for the other and 
no quantitative reliability can be guaranteed. In fact, in 
getting the estimate Xph,Fe/^ph,Sc = 0.45 we have ig- 
nored the possible differences in the p- and /- DOSs of 
the two solids. Consideration of these differences would 
lead to an inferior quantitative agreement. Also, the ra- 
tio ra{LiP) for Fe to Sc for the volumes considered is 2.35 
according to Eq.([T5]), while the value from the linear re- 
sponse calculation is 3.5. 



VII. SUMMARY 

FP-LMTO linear response calculations for the hep 
phase of Sc, based on a fixed c/a ratio of 1.592, shows 
a monotonic increase in electron-phonon coupling with 
pressure. Calculated phonon frequencies for the equilib- 
rium lattice parameter are 6-10% lower than those ob- 
tained via INS experiments.^^ The agreement can per- 
haps be improved with calculations done on a finer wave 



vector mesh. The estimated pressure, based on the lin- 
ear response results, for the instability to appear is 23- 
30 GPa. This pressure range, which can certainly be 
narrowed with additional calculations, agrees with the 
experimental observations. The F-point E2g frequency 
shows a modest increase with pressure, in qualitative 
agreement with the Raman work by Olijnyk et al^ 
Energy-optimized choice of the c/a ratio for each volume 
per atom may lead to a better agreement. 

The electron-phonon coupling constant Xph is found to 
increase steadily with pressure in the hep phase, until the 
pressure reaches a value where the hep phase becomes un- 
stable. An estimate of the electron-paramagnon coupling 
constant based on the measured temperature coefficient 
of the electronic specific heat, calculated band density of 
states and Xph suggests that the spin fluctuations at nor- 
mal pressure should be strong enough to suppress super- 
conductivity completely. At the highest pressures where 
the hep phase is still stable, the increase in Xph and a de- 
crease in Xsf in proportion to the calculated band density 
of states suggest the possibility of a very low Tc super- 
conductivity, as noted by Wittig et al^ A comparison 
of the band densities of states in the hep and the Sc-II 
phases shows that the spin fluctuation effects in the Sc-II 
phase should be reduced by a factor of two or more. It is 
argued that this suppression of spin fluctuation combined 
with electron-phonon coupling constants of a magnitude 
similar to that calculated for the high pressure hep phase 
can indeed account for the observed superconductivity in 
the Sc-II phase. 
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